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Abstract. The thermodynamic potential and thermal dependence of low lying mass spectra of 
scalars and pseudoscalars are evaluated in a generalized Nambu — Jona-Lasinio model, which 
incorporates eight-quark interactions. These are necessary to stabilize the scalar effective potential 
for the light and strange quark flavors, which would be otherwise unbounded from below. In addition 
it turns out that they are also crucial to i) lower the temperature of the chiral transition, in conformity 
with lattice calculations, ii) sharpen the temperature interval in which the crossover occurs, iii) or 
even allow for first order transitions to occur with realistic quark mass values, from certain critical 
values of the parameters. These are unprecedented results which cannot be obtained within the NJL 
approaches restricted to quartic and six-quark interactions. 
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Nambu - Jona-Lasinio (NJL) models ([1], for reviews see e.g. [2, 3]) have the very ap- 
pealing property of describing dynamical breakdown of chiral symmetry. In the present 
talk we show our recently obtained results [4], drawing particular attention to two dis- 
tinct patterns of chiral symmetry breaking (SB) and their impact on the nature, temper- 
ature values and slopes of chiral transitions. This study has its roots in the observation 
[5] that extensions of the NJL model to accomodate the approximate SU(3) flavor sym- 
metry of the u, d, 5-quarks and the Ua( 1) breaking instanton induced 't Hooft interaction 
[6], display an unstable/ metastable vacuum in stationary phase (SPA)/ mean field ap- 
proximations, and subsequent resolution of this problem by the addition of eight-quark 
interactions to the Lagrangian [7, 8]. The multi-quark interactions considered are the 
most general non-derivative chiral symmetric spin zero combinations. A set of stabiliza- 
tion conditions constrain the coupling strengths, from which the N c dependence of the 
OZI- violating eight-quark interactions is inferred. Furthermore SPA coincides then with 
the mean field approach. The characteristics of the low lying pseudoscalar and scalar 
nonets at T = have been reevaluated in the present framework [8]. One main conclu- 
sion is that identical spectra, except for the scalar singlet-octet mixing channel (strongest 
effect on the c-meson mass) can be obtained from two distinct effective potentials. They 
are generated by just changing the strengths of 4- and 8-quark (q) couplings keeping all 
other model parameters fixed. The 4q coupling regulates the curvature at the origin and 
thus determines either the Wigner-Weyl or broken phase (in absence of other interac- 
tions). Higher order interactions can however induce symmetry breaking on top of the 
Wigner-Weyl phase, i.e. a second minimum arises with a finite condensate while the 
origin keeps further its status as a minimum. This is in contrast with the case in which 



the curvature at the origin represents a maximum of the effective potential. Then higher 
order interactions will not generate further minima, but simply shift the position of the 
existing minimum. Here their action results in a perturbative effect around the broken 
phase, while in the former case they are the motor for non-perturbative SB. These start- 
ing configuartions with double vacua at T = lead to a lowering of the critical temper- 
ature. The phenomenom of multiple vacua is known to occur within several approaches 
to the QCD vacuum [9]. We show below that these patterns of SB are still present for 
realistic values of quark masses, although the origin loses of course its significance as 
a reference point for the curvature. Explicit fits reveal further that it is the 't Hooft 6q 
strength and not the 4q or 8q couplings which induce SB. This does of course not pre- 
clude the important role played by these interactions, without which neither stability of 
the vacuum nor a "twin fit" of mass spectra is possible. 

The present analysis focuses on the gap equations at finite temperature, whose solu- 
tions represent the extrema of the thermodynamic potential. The expressions were de- 
rived within a generalized heat kernel scheme [10] which takes into account quark mass 
differences in a symmetry preserving way at each order of the expansion at T = in [7] 
and we will consider from now on the isospin limit, m u = mj ^ m s 

h u + ^M u (3I -A us h) = 0, h s + -^M s (3I + 2A U sh)=0, (1) 



which must be solved selfconsistently with the stationary phase equations 



Ght + A,- + -^ hjh k + ^ hi(hj + hj + h 2 ) + ^h] = 0. (2) 



Here Ay = M 2 — M 2 , A; = M,- — m ; -, i, j, k = u,d,s (with cyclic permutations of u, d, s 
for three possible equations) and i^j^ k, Mi denote the constituent quark masses. It is 
obvious that h u = h^ for the considered case. The factors I{ are given by the average 

oo 

Ii = l - [2Ji{M 2 u ) + J { (M 2 )] , J t (M 2 ) = J-^,p(tA 2 )exp[-tM 2 l (3) 

o 
and represent one-quark-loop integrals with the Pauli-Villars regularization kernel [11] 
p(tA 2 ) = 1 — (1 + ?A 2 )exp[— tA 2 ], where A is an ultraviolet cutoff (the model is not 
renormalizable). For this case one needs only to know 

7 (M 2 )=A 2 -M 2 ln(l + ^), 7 1 (M 2 )=ln(l + ^)- x 4^. (4) 

The model parameters are the four quark coupling G ~ N^ , the 't Hooft interaction 
coupling k ~ AfT 3 , the eight-quark couplings g\,g2 (gi multiplies the OZI violating 
combination), the current quark masses m, and the cutoff A. The stability of the effective 
potential is guaranteed if the couplings fulfill the following inequality [7] 

gi>0, si+3£ 2 >0, G> — (£) , (5) 

g\ V16/ 

from which we deduce that g\ must scale at most as A^ 5 and at least as N~ 4 [8]. 
The generalization to finite temperature of these expressions occurs in the quark loop 
integrals Jq.Ji. After introducing the Matsubara frequencies [12] 



TABLE 1: Parameters of the model at T = 0. The couplings have the following units: G 
(GeV~ 2 ), k (GeV~ 5 ), g\, g2 (GeV~ 8 ), m u = nid,m s , and A are given in MeV. The values of 
constituent quark masses M u = Mj and M s are shown in MeV (only the case of global minima). 



Sets 


m u 


m s 


M u 


M, 


A 


G 


— K 


g\ 


82 


a 


5.8 


183 


348 


544 


864 


10.8 


921 


0* 


0* 


b 


5.8 


181 


345 


539 


867 


9.19 


902 


3000* 


-902 


c 


5.9 


186 


359 


544 


851 


7.03 


1001 


8000* 


-47 


d 


5.8 


181 


345 


539 


867 


5.00 


902 


10000* 


-902 



TABLE 2: The masses, weak decay constants of light pseudoscalars (in MeV), the singlet-octet 
mixing angle d p (in degrees), and the quark condensates (uu) , (ss) expressed as usual by positive 
combinations in MeV 

Sets m n rriK m-q m^ f n f% d p —(ilu)^ — (ss)i 

a 138* 494* 480 958* 92* 118* -13.6 237 191 

b 138* 494* 480 958* 92* 118* -13.6 237 192 

c 138* 494* 477 958* 92* 117* -14.0 235 187 

d 138* 494* 480 958* 92* 118* -13.6 237 192 



, oo 

J (M 2 ) ^ J (M\T) = \6n 2 T £ J** Jdsp(sA 2 )e-^+^ 2T2+ P 2+M2 \ (6) 
and using the Poisson formula 

oo oo /"-(-co 

E F(n) = £ / dxFixy 2 ™, (7) 

n—— <x> m=— oo °° 

where F(n) = exp[— s(2n + l) 2 n 2 T 2 ], one integrates over the 3-momentum p leading to 

OO I" r, 

7o(M 2 ,r) = |^p(,AV^ 2 H2E(-l)"ap(^) 



o " =1 

f2 7^^ _ _ 3T~(*sr2 T\ I2\/i2 



(8) 



Similarly one gets J\ (M 2 , T) = -dJ (M 2 , T)/dM 2 . One recovers at T = the starting 
expressions (4) and verifies also that limT^^Joi (M 2 , T) = 0. 

Using these formulae in (1), we solve the system (l)-(2) numerically, assuming that 
the model parameters G, K,gi,g2,nii,A do not depend on the temperature. As a result 
we obtain the temperature dependent solutions M,(r), representing the extrema of the 
thermodynamic potential. 

The fit of the model parameters is obtained by fixing low lying pseudoscalar and scalar 
meson characterisitics at T = 0, (stars denote input) in Tables 1-3. As already observed 
in [8], /o(600) is the main observable responsive to changes in the OZI violating eight- 
quark interaction term, diminishing with increasing strength of the g\ coupling. The 

Mi(T) are shown in Fig. 1 (sets c andd). There are either one or three (M„ =Mj,Ms), 
i = 1, 2, 3 couples of solutions at fixed values of T. For set (c) (as well as (b)) only one 
branch of solutions is physical, i.e. positive valued. The other two have negative values 
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FIGURE 1. Left: Branches of M u l '(T),Ml (T) pairs, denoting extremal points of the thermodynamic 
potential, given by the solutions of the gap equations as functions of the temperature. Solid lines (start 
at T = as deepest minima), dashed lines (start as relative minima at T = 0) and dotted lines (saddle at 
T = 0) for the parameter set (c). Only one branch is in the physical (positive mass) region (solid curves). 
Right: The same as on the left but for set (d). All branches lie in the physical region and coexist up to 
T = 7/,, from this value on only one branch survives, with much lower values of M,. 



for the light quark masses. One sees however (set c) that the onset of the transition 
occurs at a value ofT = T a for which the other unphysical two branches meet and cease 
to exist. The rapid crossover occurs in the short temperature interval 125 < T < 140 
MeV. The crossover pattern is in contrast with the St/ (3) limit case with zero current 
quark masses, where one branch collapses to the origin M u = M^ =M S = for all values 
of the remaining model parameters and T. In this case the transition is first order [13]. 

We observe however that below a certain critical value of GA (accompanied by a 
critical value of g\) one obtains, also for the case of realistic quark masses, solutions 
with all branches positive valued at any T. This is the case shown in Fig. 1 set (d). 
Two of the branches (starting from the stable minimum and the saddle solution at T=0) 
merge in the physical region at a certain Tf, and the surviving branch has a significantly 
lower mass value. This leads to the discontinuities in observables typical of first order 
transitions. The decrease in temperature observed in sets (c) and (d) is welcome in view 
of recent lattice calculations [14], obtained for finite values of the quark masses. In this 
case there is evidence that a rapid crossover occurs, as opposed to an expected first order 
transition for the massless case [15]-[18]. Lattice QCD data have not unambiguously 
settled the question about the order of the chiral transition. For physical values of the 
quark masses, calculations with staggered fermions favor a smooth crossover transition 
[16], while calculations with Wilson fermions predict the transition to be first order [19]. 
At zero chemical potential there is growing evidence that the transition is crossover, 
which would set an upper bound for the OZI- violating 8q coupling g\. 

We finally remark that the parameter set (a), without eight-quark interactions, evolves 
as function of T qualitatively as in Fig. 1, however the crossover takes place at much 
larger temperatures, T ~ 210 MeV. Also the transition is much smoother than for set (c). 

The masses of scalar and pseudoscalar mesons at finite temperature obtained for the 
set (c) are shown in Fig. 2. As can be seen there is a rapid crossover for all meson masses 
in the same temperature interval as in Fig. 1. However, neither this rapid crossover nor 
the first order transition case (d) do imply restoration of chiral or 6^(1) symmetry, 
but only the recovery of a distorted Wigner - Weyl phase, with the minimum of the 
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FIGURE 2. The masses of pion, a = ,/b (600), 77, kaon,n', ao, A^ and /n(980) from bottom to top, for 
set (c), as functions of T (all in MeV). 



thermodynamic potential shifted to finite quark mass values due to flavor breaking 
effects. 

The role played by the different multi-quark interactions can be further understood by 
analyzing the following two limits, with the parameter set (c) as starting condition. 

Case 1: We set gi = g2 = k = and remaining parameters as in (c). In this limit the 
gap equation has only one solution for the considered parameter set, thus the system is 
in a distorted Wigner - Weyl phase, 

Case 2: We set K = and all other parameters fixed as in (c). In this case there is no 
C/a(1) breaking, but OZI violating effects are present. We verify that in this limit the gap 
equation has also only one solution, being again in a distorted Wigner - Weyl phase. 

Thus the spontaneous symmetry breakdown seen in the full set (c) at T = (and also 
in sets b and d) is driven exclusively by the 't Hooft interaction strength K. We wish not 
to include case (a) in the present discussion, as it violates the stability conditions of (5). 

TABLE 3: The masses of the scalar nonet (in MeV) at T = 0, and the corresponding singlet-octet 
mixing angle G s (in degrees). 



Sets 


m a (980) 


m K*(800) 


m /o(600) 


m /o(980) 


Os 


a 


963.5 


1181 


707 


1353 


24 


b 


1024* 


1232 


605 


1378 


20 


c 


980* 


1201 


463 


1350 


24 


d 


1024* 


1232 


353 


1363 


16 



In conclusion, the present study indicates that chiral eight-quark interactions have a 
strong effect on the temperature dependence of observables described by NJL models, 
offering a plethora of solutions deeply rooted in the nature of chiral transitions. Within 
the model its origins can be traced back to the pattern of dynamical chiral symmetry 
breaking. We have discussed at length how the different patterns emerge. Mesonic 
spectra built on the spontaneously broken vacuum induced by the 't Hooft interaction 
strength, as opposed to the commonly considered case driven by the four-quark coupling, 
undergo a rapid crossover to the unbroken phase, with a slope and at a temperature 
which is regulated by the strength of the OZI violating eight-quark interactions. This 



strength can be adjusted in consonance with the four-quark coupling and leaves the 
spectra unchanged, except for the sigma meson mass, which decreases. This effect 
also explains why in the crossover region the sigma meson mass drops slightly below 
the pion mass. A first order transition behavior is also a possible solution within the 
present approach. Additional information from lattice calculations and phenomenology 
is necessary to fix finally the strength of interactions. We expect that the role of eight- 
quark interactions are of equal importance in studies involving a dense medium and 
extensions of the model with the Polyakov loop [20]. The latter is known to increase the 
transition temperature by ~ 25 MeV [14]. In relation with the NJL model the role of the 
Polyakov loop has been investigated in several papers, see e.g [21]. The present study 
can be extended likewise. In the two flavor NJL the inclusion of eight quark interactions 
has been analyzed in connection with finite T, chemical potential and Polyakov loop 
[22], where the relevance of eight-quark interactions has been reported. Eight-quark 
physics has further been explored in presence of a constant magnetic field [23], and also 
in that case it provides for a rich structure of the effective potential. 

This work has been partly supported by grants of Fundacao para a Ciencia e Tecnolo- 
gia, POCI 2010 and FEDER, POCI/FP/63 930/2005, POCI7FP/8 1926/2007. 

REFERENCES 

1. Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 124, 246 (1961); V. G. Vaks and A. I. 
Larkin, Zh. Eksp. Teor. Fiz. 40, 282 (1961) [Sov. Phys. JETP 13, 192 (1961)]. 

2. S.P Klevansky, Rev. Mod. Phys. 64, 649 (1992). 

3. T. Hatsuda, T. Kunihiro, Phys. Rep. 247, 221 (1994). 

4. A.A. Osipov, B. Hiller, J. Moreira, A.H. Blin, Phys. Lett. B 659, 270 (2008) 

5. A.A. Osipov, B. Hiller, V. Bernard, A.H. Blin, Ann. of Phys. 321, 2504 (2006). 

6. G 't Hooft, Phys. Rev. D 14, 3432 (1976); G. 't Hooft, Phys. Rev. D 18, 2199 (1978). 

7. A.A. Osipov, B. Hiller, J. da Providencia, Phys. Lett. B 634, 48 (2006). 

8. A.A. Osipov, B. Hiller, A.H. Blin, J. da Providencia, Ann. of Phys. 322, 2021 (2007). 

9. P. Bicudo, J. E. Ribeiro, A. V. Nefediev, Phys. Rev. D 65, 085026 (2002). 

10. A. A. Osipov, B. Hiller, Phys. Lett. B 515, 458 (2001); Phys. Rev. D 64, 087701 (2001); Phys. Rev. 
D 63, 094009 (2001). 

11. W. Pauli and F. Villars, Rev. Mod. Phys. 21, 434 (1949). 

12. J. I. Kapusta, "Finite-Temperature Field Theory", Cambridge: Cambridge University Press, 1989. 

13. A. A. Osipov, B. Hiller, J. Moreira, A. H. Blin, J. da Providencia, Phys. Lett. B 646, 91 (2007). 

14. Y. Aoki, G Endrodi, Z. Fodor, S. D. Katz, K. K. Szabo, Nature 443, 675 (2006); Y Aoki, G. Endrodi, 
Z. Fodor, S. D. Katz, K. K. Szabo, Phys. Lett. B 643, 46 (2006). 

15. R. D. Pisarski, F Wilczek, Phys. Rev. D 29, 338 (1984). 

16. FR. Brown et al, Phys. Rev. Lett. 65, 2491 (1990). 

17. H. Meyer-Ortmanns, Rev. Mod. Phys. 68, 473 (1996). 

18. J.T Lenaghan, D.H. Rischke, J. Schaffner-Bielich, Phys. Rev. D 62, 085008 (2000). 

19. Y Iwasaki, K. Kanaya, S. Kaya, S. Sakai, and T Yoshie, Zeit. Phys. C 71, 343 (1996). 

20. A. Polyakov, ICTP report IC/78/4 (1978), unpublished; G. 't Hooft, Nucl. Phys. B 153, 14 (1979); 
B. Svietitsky, Phys. Rep. 132, 1 (1986). 

21. K. Fukushima, Phys. Lett. B 591, 277 (2004); E. Megias, E. Ruiz Arriola, L. L. Salcedo, Phys. Rev. 
D 74, 065005 (2006); AIP 892, 444 (2007); H. Hansen, et al., Phys. Rev. D 75 065004 (2007). 

22. K. Kashiwa, H. Kouno, T Sakaguchi, M. Matsuzaki, M. Yahiro, Phys. Lett. B 647, 446 (2007); K. 
Kashiwa, M. Matsuzaki, H. Kouno, M. Yahiro, Phys. Lett. B 657, 143 (2007); Phys. Lett. B 662, 26 
(2008); arXiv:0803. 1902 [hep-ph] 

23. A. A. Osipov, B. Hiller, A. H. Blin, J. da Providencia, Phys. Lett B 650, 262 (2007); SIGMA 4, 024 
(2008). 



